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Abstract 

The Hawking radiation for massless spin-2 fields is numerically studied when the spacetime 
background is (4 + n)-dimensional Schwarzschild black hole phase. In order to check the validity 
of the Emparan-Horowitz-Myers argument, black holes radiate mainly on the brane, we assume 
that the radial equation for the massless spin-2 fields propagating on the brane obeys the master 
equation approximately. The transmission coefficient is computed explicitly by making use of the 
Hawking-Hartle theorem. It is shown that the total emission rates into the bulk are dominant 
compared to the rates on the visible brane when n > 3. However, the bulk-to-brane relative 
emissivities per degree of freedom always remain O(l) roughly. The experimental significance of 
these results in the production of mini black holes in future colliders is briefly discussed. 
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One of the most important consequence of the recent brane-world scenarios with 
large or warped extra dimensions is the emergence of low-scale (~ 1 TeV) quantum gravity. 



;he mini black holes in the future 
4J] . In this reason the absorption 



This fact opens a possibility for the copious production o 
colliders such as LHC by high-energy collision experiment 
and Hawking radiation for the higher- dimensional black holes have been extensively explored 
recently. 

Emparan, Horowitz and Myers(EHM) argued in Ref. |j] that the higher-dimensional black 
holes radiate mainly on the brane via the Hawking radiation. This argument was supported 
numerically in the case of standard model(SM) field emission by the (4+n)-dimensional non- 



rotating black holes 



0,0,3. 



EHM argument was also examined in the higher-dimensional 
rotating black hole background [9]. When black holes have angular momenta, it is well-known 
that there exist the superradiance modes. It was argued that the existence of the super- 
radiance modes may lead a different conclusion from EHM argument. However, numerical 
calculation shows that EHM argument still holds in the scalar emission by 5d rotating black 

n 

holes with two different angular momentum parameters [101]. 

In this lettter we would like to re-examine the EHM argument in Hawking radiation for 
massless spin-2 fields when the spacetime background is (4 + n)-dimensional Schwarzschild 
black hole whose metric is 



ds z = -hdt z + h~ l dr 2 + r W 



n+2 



(1) 



where h = 1 — (r///r) n+1 and the angle part <ia7 2 +2 is a spherically symmetric line element 
in a form 



dtt 2 n+2 = d6 2 n+l + sin 2 6 n+1 



d6 2 n + sin 2 6 n 



sin 2 9 2 (d9 2 + sin 2 Q x d^) 



(2) 



Since graviton is not localized on the brane unlike the SM particles, its emission spectrum 
may exhibit different behaviors from spectra for other fields. 
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The emission of the graviton into bulk was numerically explored in Ref. jll. 
Following the Regge- Wheeler method, it is well-known that the gravitational perturbations 
in the spacetime dimensions larger than four consist of three modes according to their 
tensorial behavior on the spherical section of the background metric: scalar, vector and 
teasorQ QyQ. Thus the emissl o„ spectra for the graviton carr be computed via the 
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Hawking formula 20|, |21 ] 

d 2 T BL = (n + l)(n + 4) 
duett 2 



2 n+2 7r (r i .+3)/2p 



n + 3 



-i 



O A 



EU 

A=S,V,T 



(3) 



where (3h is an inverse Hawking temperature, and S, V, and T denote the corresponding 
scalar, vector, and tensor modes. Of course, erf L is an total absorption cross section for each 
mode. 
In Ref. 



Ill ] the bulk graviton emission rate is compared to those for the SM fields propa- 
gating on the brane and concluded that the bulk graviton emissivities are highly enhanced 
with increasing n. However, the bulk-to-brane ratio for the graviton emissivities was not 
computed in the paper. Thus, the result of the paper does not lead any conclusion whether 
EHM argument is valid or not in the problem of the spin-2 field emission. In Ref. 12| the 
ratio of bulk graviton emissivity to bulk scalar was computed, which is summarized in Table 
I. Table I indicates that the emission of the graviton fields become dominant more and more 
with increasing n. In this Letter we will compute the brane decay rates for massless spin-2 
fields numerically and as a result, we will show that the bulk-to-brane ratio of the total 
emissivity becomes 0.76, 0.66, 1.59, 4.25 and 23.93 when n = 1, 2, 3, 4 and 6 respectively. 
This indicates that the total emissivities into the bulk become dominant when n > 3. Since, 
however, the bulk spin-2 fields have (n+4)(n+l)/2 polarization states while brane fields have 
only two helicities, the emission rates per degree of freedom into the bulk become roughly 
same order with those on the brane, which strongly support the EHM argument. In spite of 
0(1) roughly in the bulk-to-brane relative emissivities the dominance of the bulk emission 
rates indicates that we cannot ignore the missing energy portion in the future experiments 
relating to the brane-world black holes. 

Table I: Graviton-to-Scalar Ratio in Bulk Emission 



n 


1 


2 3 4 


6 


Graviton / Scalar 


0.052 1.48 


5.95 12.1 18.8 


34.1 



Now we would like to discuss the emission of the spin-2 fields on the brane, whose metric 
is projected from the (4 + n)-dimensional Schwarzschild spacetime (CQ): 



dsl = -h{r)dt 2 + h-\r)dr 2 + r 2 {d0 2 + sin 2 60 2 ). 



(4) 
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In fact, the spin-2 graviton is a particle living in the bulk. Thus, it seems to be ridiculous to 
consider the graviton propagating on the brane. We think this bulk-nature of the graviton is 
a main reason why the radial equations for the axial- and polar-perturbation are not uniquely 



determined 



22 



23l | when the spacetime is a projected metric (jl]). Since the purpose of this 



paper is to check the validity of the EHM argument in the spin-2 field, we will consider 
the hypothetical spin-2 field whose radial equation is assumed to be obeyed by the master 
equation as follows. 

In Ref. |24| the perturbations for the scalar(s = 0), fermion(s = 1/2) and 



vector(s = 1) fields were discussed in this background by employing Newman- Penrose for- 
malism and the following radial master equation was derived: 



A 2 Y + PA_Y - QY = 
where A± = d/dr* ± iu, A 2 = A + A_, d/dr* = h d/dr and 



dr* ^ U \h 



Q = % \A ls + (2s + n + l){ns + s + 1)(1 - h)\ 



(5) 



(6) 



with A ls = £(£ + 1) - s(s + 1). 

Although Eq.Q was derived without considering the graviton(s = 2), one can easily 



25|. In Ref. 



26j, furthermore, Eq.flSJ) is 



show Eq.(J5J) is valid for the graviton when n = 
assumed to be valid for arbitrary positive n for the graviton and derived the physically 
relevant quasinormal frequencies. However, one can show by directly applying the Newman- 
Penrose formalism that the master equation (jSJ) is not valid for the spin-2 graviton fields 
propagating on the brane. As commented above, we guess the proper radial equation for 
the graviton confined on the brane does not exist due to the bulk-nature of the graviton 
field. Since, however, the original purpose of the present paper is to check the validity of 
the EHM argument of the spin-2 field, we consider the hypothetical spin-2 field whose radial 
equation is assumed to be Eq.flS]) with s — 2. In this Letter, therefore, we will use Eq.flSJ) 
for the computation of the emission spectra for the spin-2 fields propagating on the brane 1 . 

Defining R = f~ l Y where / is defined as (1 / f)df / dr* = —P/2, one can transform Eq.flSJ) 
into the Schrodinger-like expression A 2 R = Vb r R where the effective potential Vf, r is in general 



1 Another reason for the choice of Eq.© as a radial equation is as follows. Since the quasinormal frequen- 
cies computed in Ref.[26( are physically reasonable, we can assume that Eq.© is to some extent valid 
approximately to describe the graviton propagating on the brane. In this reason Eq.([5]) can be used for 
the approximate computation of the Hawking radiation into the graviton confined on the brane. 
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complex in the following expression 

K^P + ^ + I^ + g. (7) 

Since the method for the computation of the reflection and transmission coefficients is ex- 

[ 

plained for the case of the complex potential in Ref . [251] , we will adopt the procedure of 

Chandrasekhar for the computation of the absorption and emission spectra. The solution 
convergent in the near-horizon and asymptotic regimes are respectively 

Rnh = C NH r]f s ( ) d ^ x ~ X ^ N+Pn ( 8 ) 



N=0 

00 



) x-( 7V+1 ) 



= u- s x s+1 e~ tx tn(+)X- {n+1) + C 00 u'x-' +1 e i3 ' £ t n(- 

N=0 N=0 

with x = ur, xh = urn, de t o = r (±) = 1 and 

The sign for p n is chosen by Imp n < to ensure the incoming behavior of the spin-2 wave 
in the near-horizon regime. The recursion relations for the coefficients d^N and Tjvy±) can 
be explicitly derived by inserting Eq. ((SJ) into the wave equation, i.e. A 2 R = Vb r R- 

The transmission coefficients Tbr for the complex potential (|7j) can be derived as follows. 
Firstly, we relates changes of the black hole mass(<iM) to changes of the horizon area(rfS). 
Although this relation is simply 

AG 

cffi = ^dM (10) 

for the bulk metric (pQ), where T# and G n are the Hawking temperature and (4 + Tri- 
dimensional Newton constant, it becomes slightly complicated form for our case as follows; 

d£ = - J" , dM (11) 

where fin = 1/Th, 

2 5+2n 7T 2 G n 

in ~ (n + l)"(n + 2)n n+2 (12) 
and Q n +2 is an area of a unit (n + 2)-sphere. 



Employing the Hawking-Hartle theorem 23, 123, |29| we express the variation in the area in 



terms of the variations in the spin coefficients. Finally the variations in the spin coefficients 
are identified with the perturbation in the Weyl tensor, using Ricci identities. Assuming Y = 



rh 2 R satisfies the master equation (jSJ) where the Weyl tensor \& is = R(r)S(6)e 
one can derive the transmission coefficient in the following expression 

(n + 2)a n+2 r n H - 6 iu 2 



imm—iujt 



T, 



BR 



\c 



NH 



(13) 



87rG n [uo 2 + (2n//3 H y] 

Since we adopt the unit Gq — 1, we should know the relation between G n and Go- Let 
7„ = G n /Go. Then 7„ can be numerically computed by noting that Tbr in Eq. ffTBl should 
be saturated to unity in u — > oo limit. The numerical result strongly suggests that Cnh is 
dependent on neither to nor n, which implies j n — (n + 2)Q n+2 /87r. 




Fig.1(a) 




Fig.1(b) 



FIG. 1: The w-dependence of Tbr when £ = 2 (a) and £ = 3 (b). The increasing rate of Tbr tends 
to reduce with increasing n, which indicates that the barrier heights of the real effective potentials 
become higher with increasing n. 



Fig. 1 is a plot of Tbr as a function of the energy u for £ = 2 (Fig. 1(a)) and £ = 3 
(Fig. 1(b)). As expected Tbr is saturated to unity with increasing uj. The increasing rate 
of Tbr tends to decrease with increasing n, which implies that the barrier heights of the 
real effective potentials become higher with increasing n although we do not know the exact 
expression of the real effective potential. However, the low-energy increasing rate of Tbr 
when £ = n = 2 (see Fig. 1(a)) seems to be extra-ordinarily large, which enables us to guess 
that the width of the potential barrier in this case may be narrower compared to the other 
cases. 



Once Tbr is computed, it is straightforward to compute the emission spectrum[20|, 1211 ] 

^BR 



d 2 r, 



to 3 a BR 



dudt 



7T 2 (e^ - 1) 



(14) 



where a BR is a total absorption cross section defined cjt = Yli^i^ + l)T BR /u 2 . We com- 
pute a BR numerically by making use of the quantum mechanical scattering theories with 
numerical analytic continuation, which was introduced in detail in Ref. 1CJ, l30l |. 
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Fig.2(a) m Fig.2(b) 



FIG. 2: The w-dependence of the spin-2 field emission spectra for n = 1 (Fig. 2(a)) and n = 4 
(Fig. 2(b)). The decay rates on the brane are plotted by red color. The blue line are bulk emission 
spectra for each mode. This fugure indicates that the bulk emission rates highly increase with 
increasing n compared to the brane emission rates. 



In Fig. 2 the tu-dependence of the emission spectra is plotted when n = 1 (Fig. 2(a)) 
and n = 4 (Fig. 2(b)). The decay rates on the brane are plotted by red color. For the 
comparision the bulk emission spectra for each mode are plotted together by blue color. 
Fig. 2 shows that the bulk emissivities are in general dominant in the high-energy domain 
while the brane decay is dominant in the opposite domain. This is mainly due to the power 
difference of to in the emission spectra formula defined in Eq.flSJ) and (|14p . Fig. 2 also 
indicates that the bulk decay rates are comparatively larger than the brane decay when 
n = 4. 

For the precise comparision we consider the total emission rate defined 

f°° d 2 T 

r W„ ^JJIr < 15) 

Table II: Brane versus Bulk in T tot /Tf ot 
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n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 6 


Brane 


spin-0 


1 


8.82 


37.5 


99.8 


229 


784 


spin-1 


0.45 


12.3 


65.5 


200 


463 


1600 


spin- 2 


0.052 


6.82 


81.3 


167 


244 


582 


Bulk 


spin-0 


i 


3 




99 n 




4fl8 

4UO 


spin-1 


S 


23 


3 66 


1 6 


47 9 


1 22 


824 


V 


0.23 


2.42 


12.8 


52.5 


193 


2383 


spin-2 


s 


0.026 


2.86 


27.2 


109 


312 


1976 


V 


0.026 


2.13 


21.4 


111 


438 


5117 


T 




0.22 


4.98 


45.9 


286 


6837 



The relative total emissivities for spin-0, spin-1, and spin-2 fields are summarized in Table 
II. Each total emission rate is divided by the four-dimensional scalar rate Tf ot = 2.98 x 10~ 4 . 
The abbreviations S, V and T denote the scalar, vector and tensor modes for the bulk fields 
respectively. 

Table II shows several interesting features. Firstly, the total bulk emissivities for the 
spi-2 fields become dominant when n > 3 compared to the emission rate for the spin-2 fields 
propagating on the brane. Secondly, the emission of the spin-2 fields into the bulk becomes 
dominant in the presence of the extra dimensions compared to other bulk SM fields. In the 
brane case, however, the emission rates for the spin-2 fields are not dominant. This seems 
to be due to the fact that spin-2 graviton, in general, is not confined on the brane unlike 
the SM particles. 



Table III: Bulk versus Brane per degree of freedom 



n 





1 


2 


3 


4 


6 


spin- 





1.0 


0.40 


0.24 


0.22 


0.24 


0.52 


spin- 


1 


1.0 


0.33 


0.22 


0.20 


0.23 


0.50 


spin- 


2 


1.0 


0.31 


0.15 


0.23 


0.42 


1.37 



Table III shows the spin-dependence of the bulk-to-brane emissivities per degree of free- 
dom(d.o.f.). Since the massless spin-1 and spin-2 fields have n + 2 and (n + 4)(n + l)/2 
polarization states respectively, the relative emissivities per d.o.f. can be read directly from 
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Table II. Table III indicates that the relative emissivities per d.o.f. are always less than 
unity except spin-2 case with n — 6, which supports the EHM argument. 

In this paper the emission rates for the spin-2 fields on the brane and in the bulk are 
explicitly computed. It is found that although the total bulk emissivities becomes dominant 
when n > 3, the bulk-to-brane relative emissivities per degree of freedom remains 0(1), 
which strongly supports the EHM argument. However, as indicated by Table II, the total 
missing energy arising due to the bulk emissivities is not negligible. Thus we should carefully 
consider the missing energy portion in the future experiment relating to the blane-world 
black holes. 

It is of interest to derive a real effective potentials from the master equation §5§ by 



employing the transformation theory 



2q and interpret the results of the present letter in 



terms of the potentials. It is of greatly interest also to explore the Hawking radiation for 
the graviton in the rotating black hole background which is a still open problem. We would 
like to study these issues in the future. 
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